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We derive a central limit theorem for the number of vertices of 
convex polytopes induced by stationary Poisson hyperplane processes 
in R d . This result generalizes an earlier one proved by Paroux [Adv. 
in Appl. Probab. 30 (1998) 640-656] for intersection points of motion- 
invariant Poisson line processes in R 2 . Our proof is based on Hoeffd- 
ing's decomposition of [/-statistics which seems to be more efficient 
and adequate to tackle the higher-dimensional case than the "method 
of moments" used in [Adv. in Appl. Probab. 30 (1998) 640-656] to 
treat the case d = 2. Moreover, we extend our central limit theo- 
rem in several directions. First we consider fc-flat processes induced 
by Poisson hyperplane processes in R d for < k < d — 1. Second we 
derive (asymptotic) confidence intervals for the intensities of these 
fc-flat processes and, third, we prove multivariate central limit the- 
orems for the d-dimensional joint vectors of numbers of fc-flats and 
their fc-volumes, respectively, in an increasing spherical region. 

1. Introduction. Central limit theorems (briefly CLTs) for models of 
stochastic geometry have been considered in various papers. For example, [1] 
and [24] investigate CLTs for Poisson- Voronoi and Poisson line tessellations 
in the Euclidean plane, respectively. More general CLTs for Poisson- Voronoi 
tessellations in the d-dimensional Euclidean space M. d have been established 
in [12] and [25]. In [9], normal approximations are given for some mean- value 
estimates of absolutely regular (/3-mixing) tessellations. A CLT for station- 
ary tessellations with random inner cell structures has been derived in [13]. 
Furthermore, CLTs and related asymptotic properties for the empirical vol- 
ume fraction of stationary random sets in M rf are examined in [2, 5, 16]. 



Received February 2005; revised November 2005. 
Supported by France Telecom Grant 42 36 68 97. 

AMS 2000 subject classifications. Primary 60D05; secondary 60F05, 62F12. 

Key words and phrases. Poisson hyperplane process, point process, fc-flat intersection 
process, [/-statistic, Hoeffding's decomposition, central limit theorem, confidence interval, 
long-range dependence. 

This is an electronic reprint of the original article published by the 
Institute of Mathematical Statistics in The Annals of Applied Probability, 
2006, Vol. 16, No. 2, 919-950. This reprint differs from the original in pagination 
and typographic detail. 



1 



2 



L. HEINRICH, H. SCHMIDT AND V. SCHMIDT 



A CLT for estimators of surface area densities in the Boolean model has 
been proved in [18], while in [11] CLTs for a more general class of random 
measures associated with absolutely regular germ-grain models have been 
proved. In [19] (and references therein), the reader can find a lot of further 
CLTs for empirical characteristics of Boolean models. Consistency proper- 
ties and asymptotic normality of joint estimators for the whole vector of 
specific intrinsic volumes of stationary random sets in M. d have been derived 
in [23] and [28], while uniformly best unbiased estimators for the intensity 
of stationary flat processes have been considered in [27]. 

In the present paper we prove CLTs for the number of vertices and the 
number, as well as the volume, of /c-flats {k = 1, . . . , d — 1) induced by inter- 
sections of d (resp. d — k) hyperplanes of stationary, not necessarily isotropic 
Poisson hyperplane processes in M. d ; see Sections 3 and 4. More precisely, 
we count the number of vertices lying in the <i-dimensional ball B d with ra- 
dius r > and centered at the origin o £ M. d . In addition, we both count the 
induced fc-flats hitting B d and measure their total /c-dimensional Lebesgue 
volume in B d for any k = 0,...,d—l and we study their joint behavior when 
the radius r tends to infinity. 

Noting that the number of vertices contained in B d can be expressed in 
the form of a multiple random sum running over all (i-tuples of distinct 
hyperplanes which have a common point in B d (see Chapter 6 in [17]), 
we use Hoeffding's decomposition of [/-statistics with a Poisson distributed 
number of random variables. Hence, asymptotic normality of the number of 
vertices is obtained by proving a CLT for a Poisson distributed number of 
independent and identically distributed (i.i.d.) random variables; see The- 
orem 3.1. Using a similar representation as multiple random sum for the 
number of fc-flats hitting B d and for the total fc-volume of the fc-flats in B d , 
we generalize the latter CLT in Section 4; see Theorem 4.1. Based on these 
CLTs, we obtain asymptotic confidence intervals for the intensities of the 
induced /c-flat processes and, quite naturally, are able to consider the case 
of multidimensional CLTs. 

We should mention that the normalization in our CLTs is, up to cer- 
tain constants, with respect to the d-dimensional volume of B d raised to 
the power 1 — l/(2d). We may interpret this as an expression of long-range 
dependences generated by the hyperplanes themselves. Furthermore, the 
choice of spherical sampling regions simplifies the proofs considerably, how- 
ever, most of the results remain valid for more general families of increas- 
ing convex sampling windows. If, additionally, isotropy is assumed and no 
restriction is imposed on the orientation vectors of the intersecting hyper- 
planes, this allows to determine centering and normalizing constants in the 
CLTs (i.e., intensities and asymptotic variances) explicitly. Moreover, the 
results of the present paper, together with Lemma 4.1 in [13], which states 



CLTS FOR POISSON HYPERPLANE TESSELLATIONS 



3 



that the influence of cells hitting the boundary of is asymptotically neg- 
ligible as r — > oo, it is possible to derive CLTs for ^-facets (k = 1, . . . ,d) of 
Poisson hyperplane tessellations. 

In Section 5 we reformulate Theorem 3.1 in the particular case d = 2 and 
compare this CLT with a related result obtained by Paroux [24] for planar 
Poisson line processes. Applying again Hoeffding's CLT for [/-statistics (with 
random normalization), we obtain a considerably simple proof of the CLT 
derived in [24] by the "method of moments." 

Applications for our results arise in stochastic-geometric network mod- 
eling, both in macroscopic settings like in telecommunication (see, e.g., [8]) 
and in microscopic settings like in cell biology (see, e.g., [4]). In particular, 
in Section 4.2 we show how our central limit theorems and especially our 
(asymptotic) confidence intervals can be applied in the framework of the 
so-called stochastic subscriber line model (SSLM) for telecommunication 
networks in urban environments. The SSLM is used in the context of strate- 
gic network planning and network analysis as a flexible model depending 
only on a limited number of parameters; see [7]. Figure 1 shows a realiza- 
tion of the SSLM in the case where a Poisson line process is used to model 
the underlying road system and where two types of network components 
are placed onto the lines. Besides tessellations induced by Poisson line pro- 
cesses, the class of Voronoi type tesselations is also used in the SSLM, for 
example, in order to model serving zones; see Figure 1. Therefore, we briefly 
discuss CLTs for Poisson- Voronoi tessellations in Section 6 which recently 
have been obtained in [12]; see also [25]. 

2. Preliminaries. In this section the basic notation used in the present 
paper is introduced and a brief account of some relevant notions of stochastic 
geometry is given. For a detailed discussion of the subject, the reader is 
referred to the literature, for example, [30] and [32]. Further background 
about random tessellations, flat and hyperplane processes can be found, for 
example, in [20] and [22]. 

Throughout, let [f2,cr(Q),P] be a common probability space on which all 
random objects are defined in the present paper. Let (x,y) = J2k=i x kUk 
denote the scalar product of the coordinate vectors x = (x\, . . . ,Xd) T and 
y = (yi, . . . ,yd) T in By means of the Euclidean norm || • || = \J (■, •), we 
may define the ball = {i£ W 1 : \\x\\ < r} centered at the origin and the 
unit sphere S"^ 1 = {x 6 M d : \\x\\ = 1} in M d , respectively. Furthermore, let 
S^- 1 = {(xi,...,x d ) T E S^-.Xd > 0} be the upper unit hemisphere and 
let ffc(-) denote the Lebesgue measure in M fc ; k = 0, ...,d. This measure 
will also be used instead of the fc-dimensional Hausdorff measure in M d for 
k = 0, . . . ,d — 1. As usual, z^o(-) coincides with the counting measure, that 
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Fig. 1. Realization of the stochastic subscriber line model. 



is, uq{B) = j£B. For brevity, put 

Kd = MBf) = 
where T(s) = / °° e~ y y s ~ 1 dy for s > 0. 
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d/2 



r(d/2 + i) ' 



2.1. Stationary flat processes. For each k G {0, . . . , d — 1}, let denote 
the space of all afhne /c-dimensional subspaces in R d and = {L G : o G 
L}. A (translation- invariant, resp. rotation- invariant) point process : 0, — > 
N(J- ; ) is said to be a (stationary, resp. isotropic) k-flat process in M d if the 
intensity measure A^(-) = E<I>/ C (-) of satisfies 

A fc (^\^) = 0, 

where J 7 ' denotes the space of nonempty closed sets in M. d and N^J 7 ') is the 
family of all locally finite counting measures on the Borel <r-algebra B(JF'). 
A (d — l)-flat process is called a hyperplane process. 

Provided that the intensity measure Afc(-) of a stationary /c-flat process 
is locally finite and different from the zero measure, there exists a finite 
number > (called the intensity of $^) and a probability measure 0^ on 
B{Ct) (the so-called orientation distribution of <3?fc) such that the following 
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disintegration formula 

(2.1) A k (B) = X k [ [ t B {L + x)v d - k (dx)@ k {dL) 

holds for any B G B{A k ). Here, L 1 - G C\_ k denotes the orthogonal comple- 
ment of L G C k and 1_b(-) stands for the indicator function of the set B. 
Formula (2.1) yields a simple interpretation of the intensity \ k as ratio 

. E$ fc ({l6^:LnB^0}) . 

(2.2) Afc = — — , , — for all r > 0. 

n d - k r d ~ k 

In other words, \ k K d _ k is the expected number of /c-flats hitting the unit 
ball in M. d . On the other hand, if we use (2.1) and apply Campbell's theorem 
to the stationary random measure 

&*(•)= E "*((•) n£) 

Lesu PP (<i> fe ) 

on £(M d ), where supp($ fc ) = {Ieif: $ fe (L) > 1}, we get 

(2.3) E($ k (B) = \ k v d (B) for bounded S G B(R d ). 

Hence, X k can be regarded as mean total fc-volume of all fe-flats in the unit 
cube [0,l) d . 

In the particular case of a stationary hyperplane process <I> with intensity 
A, formula (2.1) simplifies since each hyperplane H(p,v) = {x:(x,v) = p} 
can be parameterized by its signed perpendicular distance p G R from the 
origin and its orientation vector v G Sir 1 . Using the (spherical) orientation 
distribution 6 : S(S^ _1 ) -> [0, 1] defined by 

(2.4) Q d _ 1 (B) = Q({H(0,v)eC d d __ 1 :vGB}) for B G B(S ( ^T 1 ), 
formula (2.1) can be rewritten as 

A d . 1 (B) = \ f [ l B (H(p,v))d P e(dv) for B G B(A d _i). 

Alternatively, a (spherical) orientation distribution can be introduced as 
an even (symmetric) probability measure O* on B(S d ~ 1 ) which is connected 
with 9 by e*(B) = ±(0(5 n S*- 1 ) + G(-B n S^T 1 )) for B G BiS^ 1 ). The 
symmetry condition Q*(B) = Q*(—B) expresses the identification of hyper- 
planes with antipodal orientation vectors. A stationary hyperplane process 
$ is said to be isotropic (or motion-invariant) if is the uniform distribution 
on 5+ _1 or, equivalently, 0* is the uniform distribution on S^ -1 . 
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2.2. Poisson hyperplane processes. The above parameterization of hy- 
perplanes admits the representation of the (stationary) Poisson hyperplane 
process <& (with intensity A and spherical orientation distribution 0) as a sta- 
tionary and independently marked Poisson point process = Ya>i &(Pi,Vi) 
on R 1 with intensity A and mark distribution (on the mark space Si -1 ). 

Now, let k G {0, . . . , d — 1} be fixed. The intersections of d — k (distinct) 
hyperplanes Hi,..., H d _k belonging to the support of $> induce a stationary 
/c-flat process $ k which can be represented as multiple sum, that is, 

(2.5) $ k (B) = — l — J2 i B (i?in.-n^) 

[Cl >■ H u ...,H d _ k &upp(<f>) 

for B S B(A d ). If the indicator function in the latter expression is replaced 
by Vk{H\ fl ■ • ■ Pi Hd-k H (•)), we get a stationary random measure Cfc( - ) on 
B(M. d ). Intensity \ k and orientation distribution 0^ of <3?fc are given by 

(2.6) = ■ ■ ■ t B (H(0, Vl ) n • • • n H(0, v d _ k )) 

x Vd-k(vi, v d _ k )®{dvi) ■ ■ ■ @(dv d _ k ), 

for B G B(£ k ); see, for example, Chapter 4 in [30], where v d ~ k {v\, ■ ■ ■ i v d-k) 
denotes the (d — /c)-dimensional volume of the parallelotope spanned by the 
vectors Vi, . . . , v d _ k G S^" 1 . 

The intensity A^ of can be expressed in terms of the (d — k)th in- 
trinsic volume of the Steiner convex set (or zonoid) associated with $; see, 
for example, page 161 in [17]. In the isotropic case (i.e., G is the uniform 
distribution), these formulae reduce to 

(2.7) Xk=(f)^ ( ^) V fe for k = 0, 1, . . . , d - 1. 

V « J K k \dK d J 

2.3. Hoeffding 7 s decomposition of U -statistics. The proofs of the central 
limit theorems we are going to present in Sections 3 and 4 are based on Ho- 
effding's decomposition of [/-statistics which we briefly sketch subsequently. 
A more detailed discussion can be found in [15] and [31]. 

Let Xi,X2, ••• be a sequence of i.i.d. random vectors in R rf and, for any 
fixed m > 2, let / : M. md — > R be a Borel-measurable symmetric function such 
that E|/(Xi, . . . , X m )\ < oo. A [/-statistic ui m \f) of order m with kernel 
function f is then defined by 

(2.8) [/M(/) = -L £ f(X h ,...,X im ) fom>m. 

w l<i 1 <—<i m <n 
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Note that U^' '(/) is an unbiased estimator for \i = Ef(Xi, . . . , X m ). By 
elementary rearrangements, we may write Un (f) m the following form 
(Hoeffding' s decomposition): 

n 

(2.9) UtXf) ~ li = - 5>(*i) " M) + 

n r - ; 

i=i 

where 

g(x)=E(f{X 1 ,X 2 ,...,X m )\X 1 = x) = Ef(x > X 2 ,...,X m ) 

and 

4 m) a) = e (r) (i) _1 e G k ( Xil ,. . ., Xik ), 

k=2 ^ ' ^ ' \<i x <-<i k <n 
with functions G k '■ R fcrf — > R, 2 < k <m, defined by 
Gfc(xi,...,Xfc) 

= E(- 1 ) fc_j E (E/(x il ,... > xi j .,x i+1 ,...,x m )-/i). 

j=l l<ii<— <ij<k 

The crucial outcome of Hoeffding's decomposition (2.9) can be summa- 
rized in the estimate 

(2.10) E(R^\f)) 2 <^Ef 2 (X 1} ...,X m ) forn>m 

n 

and for some constant c m < oo only depending on m. The latter result, 
provided that E/ 2 (Xi, . . . ,X m ) < oo, immediately leads to Hoeffding's CLT 
for [/-statistics (see Chapter 5.2 in [31]), that is, 

(2.11) Vn~(U}r\f) ~ M) ^Af^mW^i) - /, 2 )), 

where — denotes convergence in distribution and where Af(0, m 2 (Eg 2 (X\ ) — 
fj, 2 )) is a Gaussian mean random variable with variance m 2 {Eg 2 (X{) — fi 2 ). 

3. Point process of intersection points. Let ^ = J2i>i\p l ,v 1 ) De the 
marked-point-process representation of a stationary (not necessarily isotropic) 
Poisson hyperplane process <3? with intensity A > and a nondegenerate 
spherical orientation distribution 0, that is, 0(L H S^ 1 ) < 1 for any L £ 
This assumption on ensures that each of the stationary /c-flat pro- 
cesses generated by $ has positive intensity X k for = 0, . . . , d — 1 and 
the Poisson hyperplane tessellation induced by $ consists of bounded cells; 
see Chapter 6 in [30]. 
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In this and the next section we derive CLTs for the number ^({L € 
^:InB r V 0}) of fc-flats (fe = 0, . . . , d - 1) hitting the ball B?, as well as 
for their total /c-volume contained in B® when the radius r tends to infinity. 
In the particular case k = 0, the atoms of the point process <3?o will be labeled 
by the d orientation vectors of the intersecting hyperplanes generating the in- 
tersection points. More precisely, for any r > and B S B((S'ir l ) d ), we con- 
sider the number ^o(B d x B) of those intersection points H(Pi 1 , VfJ fl • ■ - Pi 
H(Pi d ,Vi d ) in Bf for which the corresponding orientation vectors V^, . . . ,Vi d 
satisfy the condition (VL), . . . , Vua) <EB. Here, (Vfo), • • • , V(i d )) is the re- 
ordered vector (Vi 1 , . . . ,Vi d ) such that Vi n {\ ^ ■•• ^ V(n d )-> according to an 
appropriate linear ordering ■< in S^T 1 . 

The proof that ^^{Bf. x B) is asymptotically normal as r — > oo relies on 
the following basic property of the Poisson process \&. Given the number 
N r = ^([— r, r] x S^T 1 ) of hyperplanes hitting Bf, say, iV r = n, the random 
vectors = (Pi, VI), i = 1, . . . , n, are i.i.d. (and also independent of N r ) with 
independent components, where Pi is uniformly distributed on [— r, r] and 
Vi has the distribution 0. Notice that N r is Poisson distributed with mean 
2Ar, which corresponds to (2.2) for k = d— 1. 

In this way, we get that 

(3.1) ^xBjil /B(X n ,...,Xj, 

' i<ii,...,id<iV r 

where = means equality in distribution, the sum J2* runs over pairwise 
distinct indices, and 

fB((pi,v 1 ),...,(p d ,v d )} 

(3-2) 

= xiHip^v!) n • • • n F(p d , v d ) n £ r d )iB(i; (1) , . . . , v (d) ), 

where x{K) = 1 f° r K 7^ an d x(0) = 0- Since the function : (R x 
^rf-i^d _^ {0,1} is symmetric and measurable, the right-hand side of (3.1) 
divided by (^ r ) and conditioned on N r = n is a [/-statistic of order d with 
kernel function f = fs as defined in (2.8). 

3.1. Moment formulae. Since the first components Pi of the i.i.d. random 
vectors Xi = (Pi,Vi),i > 1, are uniformly distributed on [— r, r], the expec- 
tations E/b(Xi, . . . ,X d ) do not depend on r > 0. To simplify notation, we 
put 

a% d) = E(/b(Xi, . • • , X d )f B (X d _ j+1 , . . . , Xm-j)) 

(3-3) 

= %BpQ-j+l> • ■ • , X d), 
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where g B ((px,Vi), . . . , (j)j,Vj)) = Ef B ((p u vi), . . . , (pj,Vj),X j+ i, . . . ,X d ) for 

j = 1, . . . ,d. Notice that also the second moments <Xg do not depend on 
r > 0. Now we formulate a first auxiliary result. 

Lemma 3.1. For any B £ B{{S d + 1 ) d ), 

(3.4) E^ (B d xB) = ^^Ef B (X 1 ,...,X d ) 
and 

(3.5) Var*o(fl r d x B) = £ ^20^*™. 

Proof. By the symmetry of the function f B defined in (3.2) combined 
with the independence between N r and the i.i.d. sequence Xi = (Pi,Vi),i > 1, 
it is easily seen from (3.1) that 

E^ (B d X B) = ^E(N r (N r - 1) • • • (JV r - d + l))Ef B (X u . . .,X d ). 

Note that the <ith factorial moment of a Poisson distributed random variable 
is equal to the <ith power of its mean. Thus, 

(3.6) E(N r (N r - 1) • • • (N r -d+l)) = (EN r ) d = (2\r) d 

which proves (3.4). To derive a formula for the variance Vav^>o(B d x B), 
we again utilize the symmetry of f B and employ some simple combinatorial 
arguments which lead to 

>2 



E(* (^ x B)) 
1 

\d\) 2 \ 



1 e( E* f B (X il ,...,X id )\ 

\l<h,...,i d <N r / 

( 'd\ ( d 

3=0 

xE ( J2* fB(X il ,...,Xi d )f B (X id _. +1 ,...,X i2d _ j ] 

\l<ii,-,i 2 d-j<Nr 

d E(N r (N r - 1) • • • (N, - 2d + j + 1)) 

fr> m-m 2 

x E(/b(-X"i, . . .,X d )f B (X d _j + i, . . .,X 2d -j))- 
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Finally, applying (3.6) with d replaced by 2d—j for j = 0, 1, . . . ,d, and noting 
that the summand for j = in the last line coincides with (K^q(B^ x B)) 2 , 
we obtain (3.5), which completes the proof of Lemma 3.1. □ 

Notice that, as an immediate consequence of (3.5), we obtain the limiting 
relation 

(3.7) ^Va^xB) (2A)^ 

V > r~*oo r 1d-\ ((d-1)!) 2 B 

3.2. Central limit theorem for the number of intersection points. We now 
are in a position to formulate and prove a CLT for the number ^q(B^ x B) of 
marked intersection points as r — > oo, where the centering and normalizing 
constants have been derived in Lemma 3.1 and in (3.7), respectively. 

Theorem 3.1. Let B e B((S^ _1 ) d ) be chosen such that a ( ^ d) > 0. Then, 



(2Ar) d -V2V uv r ; d! 

PROOF. Note that (3.1) is equivalent to the equality ^o(B d x B) = 

( N J[)U^(fB)- For any fixed N r = n> d, the random multiple sum U^(fB) 
is a [/-statistic as defined in (2.8) of order d with kernel function fs de- 
fined in (3.2). Let n r denote the expected value MN r = 2Xr and let = 
E/b(^i, . . . , Xd). Then, Hoeffding's decomposition (2.9) yields 

*o(B? xB)- i^^E/B(Xi, ...,X d ) 

^)_^^- 1 )^( 1 - 1 )_l) w+ (^ <(/B) 

+ ( i ) f ffs(^) - , 

where gs{x) = E/^(a;, X2, . . . , X d ) for a; € [— r,r] x Since N r is inde- 

pendent of X\,X2, ■ ■ ■ , the estimate (2.10) implies 

E((i?g) 2 |iVr = n) < ^EfUXi, ...,X d ) for n > d. 
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Hence, we conclude that 

2 /„ \ 2 



E( ( N J j R$(f B ) )" = E ( d )" E ((*$)Vr = »)P(^r =n) 



n>cZ 

caE/|(Xi,...,X d ) ¥ ((N r -l 
d 2 \\d-l 

To determine the second moment of the binomial coefficient ( N JSi) , we use 
the expansions 

k 

x(x — 1) • • • (x — k + 1) = E 

i=i 

and 

where sj 1 ^ and sj 2 ^ denote the Stirling numbers of the first and second kind, 

respectively; see, for example, [26]. From (3.6), it is seen that EiV r fc is a 
polynomial of degree k in n r . Furthermore, E((iV r — l)(N r — 2) • • • (N r — d + 
l)) 2 can be expressed as a polynomial of degree 2d — 2 in n r such that 

1 v^-nV^o. 



2d-i V V d-1 



Hence, 

p 

where — ► denotes convergence in probability. Next we show that 

By virtue of — — ► 2A = and since the difference (^ r ) — can be 

written as a polynomial of degree d— 1 in N r with coefficients not depending 
on r, it is easily seen that both 

1 ( ( N r \ Ag\ p, Q 



r d-i/2 \ y d J d\ ) r - >o ° 

and 
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Similarly, after some elementary manipulations, we find that 
(N? -nf- d(N r - n^Nt 1 ) 



Nr-n r f d ^^ N ^ k _ dN d-i 



d-l/2 

N r -n. 



k=0 

/ d-2 



3/4 I r d-2 



, fc=o / 



which in combination with the previous relation proves (3.9). Combining 
(3.8), (3.9) and 

(3 10) If^lUff 
1 j r^" 1 V d- 1 J r^oo (d-1)!' 

and applying Slutsky's lemma (see, e.g., [14]), we see that the subsequent 
Lemma 3.2 completes the proof. □ 

Lemma 3.2. Under the conditions of Theorem 3.1, 

N r 

Af(P,<T%> 



/= (EsbW - ^Ef B (X u ...,X d )\ 



Proof. Let again n r = EiV r = 2Ar and let, furthermore, fJ-B = ^9b(^i) 
E/b(Xi, . . . ,Xd)- The characteristic function of 



N r 



Cr = n r 1/2 ( ^2g B (Xi) - n r fi B 

is then given by 



. i=l 



( it Nr 

Ee 1 *^ =exp(-ii^V/is)Eexp —=Y J 9B{X i 



The characteristic function on the right-hand side can be simplified by the 
fact that N r is independent of the sequence X\ , X2 , ■ ■ ■ and that the probabil- 
ity generating function Kz Nr takes the form exp(n r (z — 1)) for any complex 
z. Thus, 

Ke^ r =expfn r Efe (it/ v / ^)9s(^i) _ i_ jL gB f Xl 



or, equivalently, logEe 1 *^ is given by 

n r E( expf -^Bp!)) - 1 - -jL<?M*i) " f^M^ 

'Til- / \/W J T ZiTby 
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(ix) \x\ 

The well-known inequality | e lx — 1 — ix — \ < for any x € M combined 
with n~ 1/2 E| 5 |(Xi)| — ► gives 



logEe 1 *^ — ► ^— for t € R, 

r— >oo 2 

which is equivalent to the assertion of Lemma 3.2. □ 

4. Extensions and applications of Theorem 3.1. The results of Section 3 
can be generalized as follows. Let again = X)i>i ^(Pi,Vi) De the marked- 
point-process representation of a stationary Poisson hyperplane process <J> 
in R rf with intensity A > and with spherical orientation distribution such 
that 9(Ln S^T 1 ) < 1 for any L £ Notice that the expositions in this 

section are, for notational ease, presented for the particular case B = (S c ^T 1 ) d 
only. 

Instead of the intersection of d hyperplanes, we consider, in Section 4.1, a 
generalized version of Theorem 3.1. In particular, for any fixed k S {0, . . . , d — 
1}, we define families of statistics ^ k {B d ) and ( k (B d ) for k = 0, . . . , d — 1, 
denoting the number of fc-flats hitting B d and their total /c-volume in B d , 
respectively. These statistics, closely related to the /c-flat intersection process 
<3?fc defined in (2.5), have the form of [/-statistics of order d — k without 
normalizing factor (j_ r i.) , that is, 



\d-k 

->d 



(4.1) 



and 



(4.2) 



* k {B d r ) = <S> k ({L eA d k :LnB^ 0}) 



, /d-k > 



l<ii,...,id-fc<JVr \j = l 



Lesupp(<i> fc ) 



/d-k 



Y* Vk[f)H(X l3 )nB? 



( d l<i 1 ,...,i d _ k <N r \j = l / 

In Section 4.3 we prove a multivariate CLT for d-dimensional vectors 
consisting of these, suitably normalized, random variables. 



4.1. CLTs for point processes of k -flats. In analogy to Section 3.1, we 
first note that the expectations Ex(H(Xi) n •■• H H(X d ^ k ) n Bf) and 
r~ k Ev k (H(Xi) n • • • n H(X d -f.) n do not depend on r > since the first 
components Pj of the i.i.d. random vectors X{ = (Pi,Vi),i > 1, are uniformly 
distributed on [— r, r]. 



14 L. HEINRICH, H. SCHMIDT AND V. SCHMIDT 

As an extension of Lemma 3.1, for B = (S^T 1 ) d , we now determine the 
first-order and second-order moments of the random variables ^k(B^) and 



QkiBr) for k = 0, . . . ,d — 1. For this reason, let 

/ fd-k \ (1d,-1h-\ \ \ 



and 

/ /d-fc \ /2c2-2fe-l \ \ 

= E<£ fc (*d-k). 

where 

(4.3) &X)fc ((p, u)) = E x (#(Xi) n • • • n FCXd-jt-i) n ( P , u) n i? r d ) 

and 

(4.4) ^((p, «)) = r- fc E^.(F(x 1 ) n • • • n (X d _ fc _i) n ff(p, u) n i? r d ) 

for G [— r, r] x Notice that the second moments a^f ^ and 

a ^vk ^ °-° a ^ so no * depend on r > 0. Using this notation, we can state 
the following moment formulae. 

Lemma 4.1. For each k = 0, . . . , d — 1, 

(2\r) d ~ k 

(4.5) EV k (B?) = P(ff(*i) n • • • n F(AV fc ) n £ r d ^ 0), 

(4.6) EC fc (B r d ) = ( ?, Ar) f * Ey fc (g(Xi) n • • • n (X d _ fe ) n 



(d-k)\ 



and 



(4 7) lim Y^M^} = ^ 2d ~ 2k - 1 .Of-*) 

( 48 ) A™ r 2d-l - fc -1)1)2^ ■ 

Proof. In analogy to the proof of Lemma 3.1, we get for k 6 {0, . . . ,d- 
1} that 



E¥ fc (fl?) = e (^ r fc )e x (h(Xi) n • • • n if(x d _ fc ) n a 
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and 



E( k (B?) = E [ d _ k \ ®MH(Xi) n • • • n #(X d _ fc ) n B?). 

Applying (3.6), we obtain that 

N r \ (2Xr) d - k 



E 



d-kj (d-Jfe)!' 

which gives both (4.5) and (4.6). Furthermore, again arguing along the lines 
of the proof of Lemma 3.1, we obtain 

Var* fe (^) = £ W 



fi j\((d-k 



H2 



/ fd-k \ / 2(d~k)~j > 

e * n#(* P )ns? M D H{x q )r\B d r 

\ \p=l I \q=d-k-j+l J 



and 



VarC fe (^) = £ 



^^((d-fc-j)!) 2 

/ /d-fe \ / 2(d-k)-j \\ 

V \p=l / \g=d-fc-j+l // 

Hence, after dividing by r 2rf_2fc_1 an d r 2d_1 , respectively, and letting r — > oo, 
we get the desired relationships (4.7) and (4.8). □ 

Recall now that the random variables ^(B^) and Ck(B d ) given in (4.1) 
and (4.2), respectively, can be expressed as [/-statistics, allowing for Hoeffd- 
ing's decomposition (2.9) to be applied. Hence, we can state the following 
CLTs, the proofs of which are in complete analogy to the proof of Theo- 
rem 3.1 and are therefore omitted. 

Theorem 4.1. Let & be a Poisson hyperplane process with intensity 
A > and nondegenerate spherical orientation distribution O. Then, for k = 
0,l,...,d-l, 

4?(X) = ^I'kZ iMBt) ~ E**(2#) +m*%t h) ) 

and 

(4.10) Zfliy) = {d ~\~ l l r ~\ Ck(B d r ) - 156(1$) 



(2Ar) rf - fc -V2 
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If the stationary Poisson hyperplane process $ is also isotropic, that is, G 
is the uniform distribution on S^ -1 , the intensity X k of the A:-flat intersection 
process & k induced by $ is given by (2.7) for k = 0, . . . , d — 1. Using this in 
the following Lemma 4.2, we establish explicit expressions for the asymptotic 
variances occurring in the CLTs (4.9) and (4.10). 

Lemma 4.2. Let G be the uniform distribution on S^T 1 and let k £ 
{0,1,..., d- 1}. Then, 

(4.11) m k {Bf) = \ k K d _ k r d - k and E( k (B d ) = \ k n d r d , 

with X k given in (2.7). Moreover, 

(1|d _fc) _ (« d - fc _i(rf-fc-l)!) 2 / d\K d \ 2 / K d -l \ 2 ( rf - fc ) 



^ 4 " 12 ^ a x> k (2d-2fc-l)! U!«fc7 Vd/« d 



and 

(4-13) <7„ fc 



_ (2 fc K d _i(d-l)!) 2 / d!/c d \ 2 //c d _i 

(2d- 1)! \k\K k J [dK d 



2(d-k) 



Proof. Both mean values in (4.11) are an immediate consequence of 
(2.2) and (2.3), respectively, where only the stationarity of the Poisson hy- 
perplane process <3? is necessary. However, in case <E> is additionally isotropic, 
the intensities A^ can be explicitely determind by (2.7). To show (4.12), we 

use the relation a^f k ^ = Eg 2 k (X d _ k ), where the function g x , k ((p, v)) is de- 
fined in (4.3). Since the random vectors Afj = (Pi, Vj), i = 1, . . . , d— k, are i.i.d. 
with independent components, where Pi is uniformly distributed on [— r, r) 
and Vi has the uniform distribution G on S 1 ^ -1 , the function g x , k ((p,v)) may 
be written in the form 



9x,k((P' v )) 



1 



(2r) 



d-k-l 



r\H(p,v)nB d j dp 1 Q(dvi) 

■■■dp d _ k _ 1 e(dv d _ k _ 1 ). 

A closed expression for g Xtk ((p,v)) is obtained by an iterated application of 
Crofton's formula 

(4.14) f [ Vj(K n H(p,v)) dp&(dv) = "f- 1 {j + 1)Kj+1 V j+1 (K), 

JS 11 - 1 JR dK d Kj 
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which holds for j = 0,1, . . . ,d — 1 and any convex compact set K C M. d ; see, 
for example, Corollary 3.3.2 in [29]. Here, Vi(-) denotes the ith intrinsic 
volume for i = 0, 1, . . . , d, where Vq{K) = x{K) an d V d {K) = u d (K) for any 
convex compact K C R d . Applying (4.14) successively for j = 0, . . . , d — k — 1, 
we get 

9xMM) = ( 2r )d-fc-i ) V d . k . 1 (H(p,v)nB r ). 

Since H(p,v) n B d is a (d — l)-dimensional ball with radius ^/r 2 — p 2 , the 
invariance and homogeneity properties of V d -k-i(') yield 

V d . k ^(H(p,v) n 5 r d ) = (r 2 -^)( d - fc - 1 )/V d _ fc _i(Sf- 1 ) 

= ( r 2 _ p2)(d-fc-l)/2 d - 1 A . 

V fe / Kfc 

see [29], page 79. Summarizing the above steps, we arrive at 
(4-15) W(P ,„,) = ^_J(-^) 
for (p,u) G [—r,r] x S^T 1 and 

\k\K k 2 d ~- k ^ ) \dK d ) Jo V ^ ^ 

Finally, we get (4.12) by observing that 

(4.16) / (l-p 2 ydp=-^—{- s = 0,l,.... 

Jo (2s + 1)! 

To verify (4.13), we make use of a formula for the second moment of Ck(B d ) 
obtained in [17], page 164: 

(4 , 7) m _ E ,g^t 2 w^^ (^ 2< " H 

for fc = 0, . . . ,d — 1. From the second formula in (4.11) and (2.7), it is seen 
that the summand for j = equals (E^(i?^)) 2 and, therefore, in accordance 
with (4.8), 

(M _ fc) ((d- A: - l)!) 2 VarOCgg 

^f.fc (2A) 2rf_2fc_1 r-Sb r 2d-l 

= d!(d-l)!K 2d „ lK ^ 1 /K d -A M " 2fc " 1 

(^0 2 «2d-2«g V2d«J 
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Finally, we obtain (4.13) by taking into account the relation 

r(d)y¥ _ 2 M ~ 1 ((d-i)!) 2 

K 2d _ 2 r(d + l/2) (2d-l)! ' 

which follows from Legendre's duplication formula 

(4.18) 2 2 *- 1 r( a + |)r(a) = V5Fr(2«) 

applied to the integer s = d. This completes the proof of Lemma 4.2. □ 
Notice that (4.5), (4.6) and (4.11) yield simple relationships between 

P(H(Xi) n • • • n H(x d _ k ) n s r d + 0) , E^ fc (iT(Xi ) n • • • n if (x d _ fc ) n £ r d ) , and 

X k for fc = 0,l,...,d — 1. Furthermore, in case <& is additionally isotropic, 
one can use (2.7) to get 

F(H(x 1 ) n • • • n H(x d . k ) n Bf. + 0) = 2k ~ ddl KdKd ~ k (Kd ' ' 



k\ K k V dn d 

r)k—dj\ 2 / \d—k 

Ev k {H{X{) n • • • n H(x d . k ) n B d T ) = r k . 

k\ K k V dK d J 

Notice that these formulae comply with results in [17], pages 160 and 161. 
Also, if we replace in the proof of (4.12) the function g Xtk {-) by g Ujk {-) defined 
in (4.4) (with V k = v k instead of Vq = x), it easily seen that 

(4.19) ^(p^))^^——^—^ (1--) 

for (p,v) £ [—r,r] x S"^ -1 , which confirms once more (4.13) without using 
Matheron's formula (4.17). On the other hand, regarding (4.17) for k = 
as a sum of power functions in Ar, we are able to determine the pair 
correlation function go(r) of the stationary and isotropic point process 
as a polynomial of degree d — 1 in (Ar) -1 ; see also [12]. More precisely, 
putting go(r) = 1 + J2jZi c dj(^ r )~^ and utilizing the relationship 

r2r 

Y&i^ Q {Bf.) = dK d \l / ^(S r d n(^ + (n,0,...,0)))(« ? o(w)-l)n d-1 du 
Jo 

+ A K d r d 

(see [32], page 131 for details), we get c dj = (^(^^(-gfL)* by com P ar " 
ison of coefficients. Here, we used (2.7) for k = together with 

u d (Bf n (Bf + (u, 0, . . . , 0)))u d - j ^ du 



d-j-l, 1 f f dxd v 



dn d Jb{Jb{ H^-ylP 
(d-j){d-j + 1)/^-, 

for j = 1, . . . ,d — 1; see [29], page 177. 
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4.2. Asymptotic confidence intervals for the k-flat intensities. Let $ be 
a stationary Poisson hyperplane process with intensity A > and nondegen- 
erate orientation distribution G. From the view point of spatial statistics, 
Afc,r = ^k{Br)/vd-k{Bfr k ), as well as X k>r = C, k {Bf) / v d {B%) , are unbiased es- 
timators for the intensity X k of the stationary fc-flat intersection process gen- 
erated by $; see Section 2.1. We mention that both estimators are strongly 
consistent since $ is ergodic and even mixing; see [30], Chapter 6.4. 

If <I> is additionally isotropic, then we know from (2.7) that X k = ad,kX d ~ k . 
Together with (4.7) and (4.12), as well as (4.8) and (4.13), we obtain 

km r l &TX X = * 2d - 2k - l 4 kb 2 d - k and km [J^^l = A M_2fc_1 a^ k b\ 

r->oo 4(d - k) 2 d,k d k r->oo 4(d _ jfc)2 d > kd 

for k = 0, . . . , d — 1, wkere 

ad,k= , — -j — and = — / - for j = 1, 2, . . . . 

\kjK k \dK d J y/2{2j - ly.Kj 

By means of (4.18), one can verify tke inequality bj < bj+i for j = 1, 2, . . . , 
wkick in turn implies tkat 

km r Var X k r < km r Var X k r for fc = 1, . . . , d — 1. 

j — >oo ' r— >oo 

Tkerefore, we prefer the estimators A^ r to construct confidence intervals 
for Afc. Notice that efficiency and other optimality properties of intensity 
estimators for stationary fe-flat processes observed in fixed convex sampling 
windows have been studied in [27]. 

By Theorem 4.1, the estimators X k)J . and X kt r are asymptotically normally 
distributed. For example, togetker witk tke above abbreviations, (4.9) can 
be formulated as follows: 

~ Xk) ^ AA(0,4(d - kfxt^alf-^U). 

Next, we apply a variance-stabilizing transformation f(x) for x > to tke 
latter CLT suck tkat y/r(f(Xk, r ) ~ /(Afc)) kas a Gaussian limit witk mean 
and variance 1; see, for example, [3] for details. It is easily ckecked tkat 
f(x) = (x j 'ad l k) 1 ' 2 ^ d ~ k ' /bd-k is a suitable ckoice for suck transformation, 
wkick gives rise to a 100(1 — a)% (asymptotically exact) confidence interval 
45(a) for A fc . 

Corollary 4.2. Under the above assumptions, 

^ a -^ d - k) \x]/^ k)) - X l '^ d - k)) ) M(0,1), 
bd-k ' ' r ^°° 
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which is equivalent to P(Afc G ljfj.(a)) — ► 1 — a for any a G (0, 1) and k 
0, . . . ,d — 1, where the interval ljfl(a>) is given by 

2(d~k) 



l/(2(d-fc)). 



z l-ct/2 



(Afc, r ) 



l/(2(d-fc)) 



a 



+ 



l/(2(d-fc)) 



N 2(«i-fc)-, 

2: l-a/2 



anc? ^i_ Q /2 denotes the (1 — a/2)-quantile of the standard normal distribu- 
tion. Moreover, relation A = (Xk/o-d^) 1 ^^^ permits to transform the inter- 
val ljfl(ct) into the interval 



4%) 



"~7T Zl -«/2 . 



\kA 1/i2(d - k)) , bd~k 



®d,k 



+ 



- z l-af. 



2i 



which covers the intensity A of the hyperplane process $ with probability 
1 — a as r — > oo . 



The results of Corollary 4.2 can be applied, for example, in the following 
context of telecommunication network modeling. 

Figure 2 shows the road system of Paris, where the real data, given by 
(x, y)-coordinates describing the location of points which are interpreted as 
intersections of roads, have been connected in order to form (nonconvex) 
cells. Beyond, each data point is equipped with a mark displaying the type 
of road to which it belongs, for example, main roads and side streets. The 
main road system is of prior interest in the civil engineering part of strategic 
network planning since these roads gather extensions of cables from different 
side streets. Therefore, expensive network components and technical devices 
are placed along main roads and capacity planning and cost analysis for the 
cables running along these types of roads is an important task. 

The SSLM introduced in Section 1 allows us to model telecommunication 
networks. An important part of the SSLM is the geometry model, that is, the 
model that represents the infrastructure or road system of a certain (urban) 
environment. In [8], a fitting procedure, based on minimization of distance 
measures and Monte Carlo test techniques, is introduced, where tessella- 
tion models are fitted to real infrastructure data. It turns out that urban 
main road systems are often best represented by a Poisson line tessellation 
(Poisson line process). 
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(a) (b) 

Fig. 2. The road system of Paris: (a) city of Paris; (b) extracted data. 

From the civil engineering point of view, it is desirable to get information 
about the network's structure without having to apply too sophisticated 
models or methods. Furthermore, it is preferable for telecommunication en- 
gineers to get such information in the form of worst case and best case 
scenario values or to test their predictions. Our (asymptotic) confidence in- 
tervals, and, based upon them, our (asymptotic) tests, contribute to such a 
risk analysis of structural parameters of the network. 

To take an example based on the observation of pointwise real data in 
a sufficiently large (spherical) sampling window B% and after application 
of the fitting procedure, it is possible to provide lower and upper bounds — 
b~(a) and 6+ (a), respectively — for the mean total length A of roads per unit 

area. To this end, consider the confidence interval Jg 2 r (a) = [b~ (a) , 6+ (a)] 
(from Corollary 4.2 for d = 2 and k = 0), where 




This means that, based only on the knowledge of the number of road 
crossings within a large (spherical) region, b~ (a)v2(W) and b+(a)v2(W) 
provide lower and upper bounds for the mean total length of the main road 
system intersecting a certain sub region W. These bounds can then be used 
to estimate lengths of cables and, in subsequent steps, to determine costs, 
as well as capacities of connection. 

Moreover, based on the above confidence interval, one is able to test the 
null hypothesis Hq : A = A* , where A* denotes some specified value of the 
mean total length of lines per unit area, versus the alternative hypothesis 
Hi : A A*. In the context of telecommunication, A* can be interpreted as the 
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ratio of the total length of cables and the area (of an urban district) in which 
the cables are laid. This value for the cable length is often determined by 
practitioners according to their own rules of thumb. Based on the confidence 

(2) 

interval Jq t .(a) (r large enough), where a is a suitable significance level, Hq 
would be rejected if 

^ 1 / 2 \ 4 ^ 1 / 2 x 4 

A , r < - I \/A* - Zi_ Q / 2 -== ) or A 0j r>-(v'A* : + ^l-a/2- 



vr V 7rv3r/ vr \ 7rv3r 

4.3. Multivariate CLTs. In this section we extend the results of Sec- 
tion 4.1 by establishing multivariate CLTs which describe the joint asymp- 
totic behavior (as r — ► oo) of the closely correlated random variables Zq^(x)i 

• • • , zf-iM), as well as of Z$(v), Z d \ r (v), defined in (4.9) and (4.10), 
respectively. To begin with, we define the mixed second moments 



(M-M-0 =EJx[ P) H(X t ) n 5? M fl H W n Bf. 



\ \i=l / \ i=d—k 

= ^(g x A x d-k)g x ,i( x d-k)) 

and 

/ /d-k \ /2d-k-l-l 

a u,ki k,d ~ l) = r~ {k+l) E ( ffe ( H H(Xi) nBrjvA fl H(Xi)nBl 



n i i o d 

i=l / \ i=d—k 



= ^(gu,k( x d-k)gu,i(x d -k)) 

for k, I = 0, . . . ,d — 1, where the functions g Xj k(-) and g v ,i{-) are defined by 
(4.3) and (4.4), respectively. The following result generalizes Lemma 4.1. 

Lemma 4.3. For any k,l = 0, . . . ,d — 1, 

(4.20) Hm E(Z« (x)<? (x)) = M ~° 
and 

(4.21) lim E(4>)^(*/)) = a^ 4 ^- 

Proof. Both relations (4.20) and (4.21) can be shown by similar com- 
binatorial arguments as employed in the proof of Lemma 3.1 to verify (3.5). 
Using the definition of zjfj.(x) as given in (4.9), we easily get that 

¥(7^MZ^M) - {d-k-l)\{d-l-l)\ 
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According to (4.1), we may write for 0<k<l<d— 1 that 

E(* fc (5 r rf )^(a?)) 

_ E(N r (N r - 1) • • • (N r - (2d - k - I - j) + 1)) 
~~ j\(d-k-j)\(d-l-j)\ 

/ (d-k \ / Id-k-l-j \\ 

xE X D ff W n5 ? X fl H{X q )^Bf . 

V \p=l / \g=d-fc-j+l / / 

Using (3.6) and the fact that the summand for j = equals E 1 I>/ C (£^?)E , I';(I^), 
we arrive at 

E(z#(x)4?eo) 

(d-fc-l)!(d-i-l)! 



x 



E 



(2Ar) M - fc - / " 1 
(2Ar) i 



d-i /n\ \2d-k-l-j 



<j\(d-k-j)\(d-l-j)\ 

/ /d-k \ / 2d-k-l-j \\ 

xE X fl^nflf x D H{X q )KB d r . 

V \p=l / V g=d-fc-j+l / / 

Hence, dividing by r 2d -k-i~i an j letting r — ► oo immediately yields (4.20). 
The proof of (4.21) is completely analogous. □ 

Theorem 4.3. Let the assumptions of Theorem 4.1 be satisfied. Then, 

(4-22) (4%))tl r ^(oMx)) 

and 

(4.23) {Z $(v))tl r -±jV(oMv)), 

where M(o,T,(x)) and A/"(o, £(z/)) are d-dimensional Gaussian mean o 

(i,d-k,d-i)\d- 

T X,k,l >k,l=0 



(0, . . . , 0) T random vectors with covariance matrices S(x) = (o - ^ 1 ^ ^)^ ^ 



and S(i^) = (cr^'hi k ' d )f u>ii/i entries given by the limits (4.20) and (4.21), 
respectively. 

Proof. Recall that due to the well-known Cramer-Wold device, the 
multivariate CLT (4.22) is equivalent to the one-dimensional CLT 

(4.24) £ t k Z$ ( X ) ^ Af(Q, t T Z(x)t) , 

k=0 
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for any t = (i , . . .,t d ~i) T €R d \ {o}. This means that the proof of (4.22) 
can be put down to the case of the (one-dimensional) CLTs considered in 
Theorems 3.1 and 4.1. First, using (2.8), (4.1) and (4.9), we may rewrite the 

linear combination J2k=o^kZjfl(x) as follows: 
V^ 1 * y(d) ( ^v* (rf-fc-1)! 

l^^^rKX) - 2^ ( * (2Xr)d - k -l/2 
k=0 k=0 ^ ' 

xK"w(A)-f$f^) 

Next, we apply Hoeffding's decomposition (2.9) to the random [/-statistic 

Ujy r k ^ (x) with kernel function x(Hj=i H(Xj)n Bf.) and proceed in the same 
manner as in the proof of Theorem 3.1. In view of the limiting relations 
(3.8) and (3.9) with d replaced by d — k for k = 0, 1, . . . , d — 1 and combined 

with Slutsky's lemma, we recognize that the weak limit of J^t^o^^^ifr (x) 
coincides with that of 

g) 4t)^-V3 ( (I WW) - 

as r — > oo. Finally, by means of (3.10) with d again replaced by d — k for 
k = 0, 1, . . . , d — 1, it remains to show that 

1 d-l I N r \ 



d-l 



±+N 0, £ tfctjE^fcCXi)^,,^)) . 



k,l=0 



However, the latter CLT is obtained by proving Lemma 3.2 once more for the 
function Y^t=o^ k 9x-k(') instead of <?b(-)- To show the second assertion (4.23), 
we only need to repeat the just finished proof with the kernel function 



-^kin-ZiHix^nB?) 

proof of Theorem 4.3. □ 



' : Vk(r\j=i H ( x j) n B r) and 9v,l{-) instead of g x ,i(-)- This completes the 



As in Lemma 4.2, the additional assumption of isotropy allows to com- 
ite explici 

(l,d-k,d-l) 



pute explicit formulae for the mixed second-order moments cr^ti k,d ^ and 



Lemma 4.4. Let O be the uniform distribution on 1 and let k,l G 
{0, . . . , d — 1}. Then, 

{l,d-k,d-l) _ (<flK d ) 2 Krf_fc_iKrf_;_i K 2 d-k-l-l ( X 2 



cr 
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(4.25) 

' ( M -oB((2rf-A;-Z)/2,(2d-fc-0/2) 

where ~B(s,t) = /q x s ~ 1 (l — x) t ~ 1 dx = r(s)r(t)/r(s+i) denotes Euler^s Beta 
function, and 

2d-k-l 

■k,d-l) _ \^d^d-l u -\ u ' ~ ± J'-J * I ^d-l 1 
~ kW.K k Kl(2d - 1)1 \dK d 

(4.26) 



a,d-k) n,d-i) 



Proof. Both (4.25) and (4.26) can be obtained using the shape of the 
functions Qx,k{') and 9v,k{') derived in the proof of Lemma 4.2. By (4.15) 
and the above definition of tf^jjf j k ' d ^ , we get that 

a (i,d-k,d-l) = — J^J^_ i g Xtk ((p, v ))g Xt i{(p,v))e(dv)dp 

d\K d Kd-k-l dlK d K d -l-l (^d- " 



k\K k 2 d ~ k ~ x Uk^- 1 ' 1 V dK d 

x [\l-pZ)W- k - l -Wd P . 
Jo 

Thus, by noting that 

[\l _ p2)(2d-fc-I-2)/2 dp = «2d-fc-t-l 

Jo 



2K2d-k-l-2 



— k—l — 2g 



2d-k-l 2d-k-l 



(see, e.g., [29], page 80), we obtain the first part of (4.25), where the second 
identity in the previous line turns out to be a simple consequence of (4.18) 
for s = (2d — k — I)/ 2 and the very definition of the Beta function. The 

second part of (4.25) is seen by inserting the variances a^ k d k ^ given by 
(4.12) combined with B(d - k,d - k) = ((d - k - l)\f/(2d- 2k - 1)! for 
k = 0, 1, . . . , d — 1. Likewise, using (4.19), we get that 

CT S5~ M_/) = tt f I ,,9v,k((p,v))g v>l ((p,v))Q(dv)dp 



d\K d K d _i d\K d K d _i ( K d _ 



\ 2d-k-l f l 

2? f^-^ 1 *- 



k\K k 2 d - k - 1 l\Ki2 d ~ 1 - 1 \ dn d 

Hence, taking (4.16) for s = d — 1, the first part of (4.26) is shown and the 
second equality is immediately seen from Lemma 4.2. □ 
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Corollary 4.4. The covariance matrix X(x) possesses always full rank 
d, whereas the rank of the covariance matrix equals 1 for any dimension 
d>l. Moreover, 

4.27 . . fc ' r /orO<Jfe<Z<d-l. 



Proof. Notice that X(x) possesses full rank if this matrix is strictly 
positive. This, however, can be seen since 

d-1 



k,l=0 



2d-k-l 1d-k-l 




x)) 



(d-k- 




for any (to, . . . , £<2-i) € R \ {o}, which means that the symmetric matrix 
with entries B((2d -k-l) /2, (2d - k - I)/ 2) for k, I = 0, . . . , d - 1 is strictly 
positive. In order to show that the rank of is 1 for any d > 1, we 

only need to observe that the second equality in (4.26) implies that each 
entry of the asymptotic covariance matrix of the normalized random vec- 
tor (zjfl (v) I J a^jf ^ )t=o e q ua l s 1- Then, (4.27) is a consequence of the 
structure of S(i/). □ 

The somewhat surprising result (4.27) states that the variance of the 

difference of any two components of (zjfj. (v) / \J a^f k ^ )^~g tends to zero 
as r — > oo. Together with Slutsky's lemma this allows for the conclusion that 
the normal convergence in (4.10) for a single component, ZqI(v), say, implies 
asymptotic normality of the other components. Thus, relation (4.27) can 
be interpreted as a kind of asymptotic second-order relationship for /c-flat 
processes induced by stationary and isotropic Poisson hyperplane process. 
It would be of interest to see whether there is a pure geometric reasoning 
for (4.27). 

5. A special review of the planar case. Throughout this section we as- 
sume that d = 2. Assuming isotropy of the underlying stationary Poisson line 
process, we first present a short review of Theorem 3.1 for the case where 
the ordered angles of the orientation vectors of intersecting pairs of lines 
are situated within a certain rectangle. In a second part of the present sec- 
tion we look at another type of a CLT for Poisson line processes, proven by 
Paroux [24], where the normalization is random. Applying directly Hoeffd- 
ing's CLT (2.11) for L^-statistics, we provide a new proof of Paroux's CLT 



CLTS FOR POISSON HYPERPLANE TESSELLATIONS 



27 



with random normalization, which has been derived in [24] by the "method 
of moments." 

5.1. Planar moment formulae. We consider the marked-point-process 
representation ^ = J2i>i \p 1 ,V 1 ) °f a planar stationary and isotropic Poisson 
line process $ with intensity A. In this special case each orientation vector 
Vi £ S\ is completely determined by the angle Tj between the unit vector 
Vi and the x-axis measured in anti-clockwise direction. Owing to isotropy, 
the angles ri,r 2) ... are independent and uniformly distributed on [0, 7r]. 
Therefore, supp( ( I > ) consists of parameterized lines ^(p^r 1 ;) i n ^ 2 defined by 
t(Pi,Ti) ={(x,y) £R 2 :xcosTi + ysmTi = Pi}. 

For r > fixed, ^o(B 2 x B(a,b)) is the random number of those inter- 
section points ^(p ni ,r ni ) ni {Pn 2 ,r„ 2 ) in B r for which (r (m) ,r (n2) ) G B(a,b), 
where B(a,b) = [0,a] x [0,6], a < b, is a rectangular subset of [0,7r] 2 . Recall 
that V (B? x B(a,b)) has the same distribution as the random double sum 

(5.i) E /Bc^C^.ro.CPj,^)), 

l<i<j<N r 

where, as in Section 3, the random vectors (Pi, Ti), (P 2 , ^2), ■ ■ ■ '■ - > r] x 
[0,7r] are independent and uniformly distributed on [— r, r] x [0,7r] and where 
N r : f2 — > {0, 1, . . .} is a Poisson distributed random variable with expectation 
2Ar independent of the (p,rj)'s. The function fB( a ,b) given by 

Pi>7i) ' f -(P2,72) 

is symmetric since (7m , 7( 2 ) ) are lexicographically ordered, that is, (7m , 7( 2 ) ) 
(71 A 72, 71 V 72). With the abbreviations 

/x B(0i6) =E/ B(ai6) ((p 1 ,r 1 ),(p 2 ,r 2 )) 

(5.2) 

= p^pia) n ^ ( p 2i r 2 ) n J3? ^ 0, < Ti a r 2 < a, ri v r 2 < b) 

and o^^^ defined by (3.3), Theorem 3.1 claims that 

(5 - 3) jdw{^ {Bl x B{aM) " r ^^°AZ^ 



where 



1 f a f b 

HB(a,b) = 4^J J \ Sill ( n ~ V ^ dudv 



(5.4) 

= — (2a — sina — shift + sin(6 — a)) 

47T 
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and 

2 

a B{a,b) = ^3 ( 5 ° + 7 ( sin ( fo - s ' mb ) 

(5.5) + (b + s'mb) cos 6 

— (b — a + sin(6 — a)) cos(6 — a)). 

To verify (5.4), one can apply the general mean value formula (2.6) to- 
gether with (2.7) for d = 2 and k = 0. However, we use (3.4) and a more 
direct approach to calculate the probability [J>B(a,b)- By definition (5.2), we 
obtain 

^ f"K rn rl rl 
VB(a,b) = ^2j o J Q j /_ ij /_ i l{x2+^<l}(Pl,7l,P2,72) 

x l[o,o](7i A 72) 

x l[ 0j6 ] (71 V 72) dpi dp 2 cfryi d^ 2 , 

where (x,y) denotes the intersection point of the lines ^( Pli7l ) and ^(p 2)72 ), 
that is, 

—pi sin 72 +P2 sin 71 



y 



sin(72 - 71) 
-pi cos 72 + P2 cos 71 
sin (72 - 71) 



Thus, by some elementary manipulations with trigonometric functions, we 
arrive at 



1 

27 2 7 



a r b rl 



J-l 



yl sin (72 - 7i ) I <^7l ^72 , 



which combined with ^/ 1 — p 2 dp = ^ confirms (5.4). To determine the 
asymptotic variance c^Ub) a PP earm g m (5-3), we argue similarly. By (3.3), 

^ ) = E (/w)(( J Pi,r 1 ),(P2,r2))/ W )((P2,r2),(P3,r 3 ))) 

^ rn f"K r7V rl rl rl 

= m Jo Jo J_J_J^+y^M^P^) 

Xl{ x 2+y-2<l}{P2, 72,P3,73) 

x%o](71 a 72 A 73) 

x l[ 0j6 ] (7! V 72 V 73) dpi dp 2 dp 3 dji d-y 2 d-y 3 
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4 pb pb pb pi 

= 8^3 Jo Jo Jo J i^ ~ Sin< ^ 72 ~ 7l ^" Sin ^ 72 ~ 73 ^' dP2 
x l[o,o] (7i A 72 A 73) efyi CZ72 ^73. 

Using that J_ 1 (l — p 2 ) dp = I, a somewhat lengthy computation of the three- 
fold integrals with respect to 71,72,73 in the latter line leads to formula (5.5). 

5.2. A CLT with random normalization. The two-dimensional version 
(5.3) of our CLT for intersection points of Poisson line processes in M 2 has 
close connections to a CLT with random normalization; see Theorem 3.1.3 
in [24]. The proof given in [24] is based on the well-known "method of 
moments." The following Theorem 5.1 states the assertion of this CLT. 
However, using Hoeffding's CLT (2.11) for [/-statistics, we obtain a much 
shorter proof. 



Theorem 5.1. For arbitrary fixed a, be [0,7r] with a <b, let 

( 5 - 6 ) Z B(a,b) = (J^V, - l))3/4 ^2 ^ (fB(a,b)(Xi,Xj) - (I B (a,b))- 
Then, 

_!> \f(o J 1 ' 2 } ,,2 \ 
Z B(a,b) r ~Z^ JV ( U > a B(a,b) ~ A*B(a,6)Js 

(1 2) 

where f^B(a,b) an d a B(ab) are 9^ ven by (5-4) and (5.5). 

(r) 

Proof. To begin with, we rewrite Z B l b -, as 

(r) Nr /2 ( 1 \ 1/4 / 2 \ 

z B ( a,b) = ^r{ 1 -w r ) \N r (N r -i) ^^y*™ {Xu Xj) " ma,6)) J 

where f B {a,b)(%i,%2) is a measurable, symmetric function in x\,X2 £ [— r,r] x 
S+ satisfying E/ B ( a)6) (Xi,X 2 ) = (J, B (a,b)- Applying Hoeffding's CLT (2.11) 
for [/-statistics yields 

(5-7) ^uW(f Bia>b) ) -^(0, Varte^CXx)), 

where the [/-statistic f/r 2 , (f B ( a ,b)) is given by 

^[„ P ] (/b(o,6) ) = r n/r j _ ^ ^ (/b(o,6) ( *t . *i) ~ MB(o,6) ) » 
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with gB( a ,b)( x ) — B(a,b){ x 1X2) i where [n r ] stands for the integer part of 
n r = ¥.N r = 2Xr. Taking into account (3.3) for d = 2,j = l and (2.9), we get 

Var 5s(0)6) (Xi) = o£g 6) - fi% M . 

p p 
Since N r — ► 00 such that N r /\n r ] — ► 1, Theorem VIII. 7.1 in [261 tells 

r— *oo r— >oo 

us that in (5.7) the deterministic integer [n r ] can be replaced by the ran- 
dom number N r without changing the limit M(0,Vai g B ^ ab ^(Xi)). Finally, a 
straightforward application of Slutsky's lemma completes the proof of The- 
orem 5.1. □ 

Notice that the asymptotic variance cr^^ft) ~~ ^B[a 6) °^ * ne number of in- 
tersection points under random centering in (5.6) is always smaller than the 
asymptotic variance o B ) ^ obtained in (5.3) under deterministic centering. 

Furthermore, instead of imposing conditions on the angles of the orienta- 
tion vectors of the intersecting lines £( Puli ) and £( pj ^^, in [24] the two angles 
at the intersecting points Xj and Xj of ^( Pii7i ) and £^ pj ^^, respectively, with 
the x-axis are considered. More precisely, the ordered pair of angles (a,/3) 
is considered, where a denotes the angle at Xj A xj between the x-axis and 
the half-line from ^( Pii7i ) C]£^ Pj ^.^ to the intersection point with the x-axis 
and, in analogy, (3 is the angle at Xj V xj. 

6. CLTs for Poisson Voronoi tessellations. In this section we give a brief 
overview of CLTs for Poisson- Voronoi tessellations (PVTs for short) in M. d , 
that is, we consider a stationary PVT generated by a stationary Poisson 
point process ^ in M. d with intensity A. Notice that each cell of the Poisson- 
Voronoi tessellation $ is defined as the closure of the set of all points in M. d 
which are closest to the nucleus of this cell, being a point of 1 I / . Moreover, 
it can be shown that each cell is a convex (d-dimensional) polytope. 

Let ^0 denote the stationary (and isotropic) point process of vertices 
of the cells induced by $ and let Ao denote its intensity, that is, Ao = 
Efo([0,l] d ). It is well known (see, e.g., [21]) that 

2 d 7T d ~ 1 K d 2 ( Kd-1 

(6.1) A = c d X with c d 



d+l K d 2_ 1 V K-d 

Consider a convex averaging sequence W n of sets in R d , that is, the sets 
W n are compact convex, increasing and contain a ball with unboundedly 
growing radius; see, for example, [6]. Since the (5- mixing coefficient of a 
stationary PVT is exponentially decaying (see [9]), a CLT for the random 
number of vertices ^o(W n ) of a stationary PVT within W n can be derived 
as n — > 00. More precisely, in [12] it is shown that 

(6.2) VdiWu)- 1 ' 2 ^^) - \ u d (W n )) -±> Af(0, A (l + c d a 2 d )), 
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where a\ is a constant only depending on the dimension d which is express- 
ible in terms of multiple integrals. A detailed discussion including numerical 
computations for the cases of d = 2 and d = 3 can be found in [10] and [12], 
respectively. The rounded values obtained there are o\ = 0.5 and o\ = 5.084, 
together with C2 = 2 and C3 = 6.768 from (6.1). 

Notice that Ao, n = ^> o{W n ) / Vd(W n ) is an unbiased estimator for the in- 
tensity Ao and, by (6.2), An jn is asymptotically normally distributed with 

mean Ao and asymptotic variance Ao(l + Crf<r^). A simple transformation 

— p 

using (6.2) and An n — ► An yields 

' n— »oo 

Similar to the confidence intervals derived in Corollary 4.2 for Poisson hy- 
perplane processes, the CLT in (6.3) enables us to construct a 100(1 — a)% 
(asymptotically exact) confidence interval for A. Indeed, for any a G (0,1), 
we have (for large enough W n ) that A 6 I n ( a ) with probability 1 — a, where 
I n (a)c d v d {W n ) is given by 

'(/w- ^aA+^I) 2 ' (a/wo + ^aA+^I) 2 " • 

Notice that an alternative and general approach to CLTs in the Poisson- 
Voronoi case is given in [25]. Among other results, a CLT for a quite large 
class of functionals related to Poisson- Voronoi tessellations is obtained there. 
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